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Abstract 

We report the results of a new calculation of soft-gluon corrections in B — > X s £ + £~ 
decays. In particular, we present the first calculation of bremsstrahlung and corre- 
sponding virtual terms to the lepton forward-backward asymmetry, which allows 
us to systematically include all contributions to this observable beyond the lowest 
non-trivial order. The new terms are important, for instance the position of the 
zero of the asymmetry receives corrections of 0(10%). Using a different method, we 
also provide an independent check of recently published results on bremsstrahlung 
and infrared virtual corrections to the dilepton-invariant mass distribution. 
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1 Introduction 



Flavour-changing neutral-current (FCNC) processes are a very useful tool to understand the 
nature of physics beyond the Standard Model (SM). The stringent bounds obtained from 
B — > X s 7 on various non-standard scenarios (see e.g. |], |2|, [| f|) are a clear example of the 
importance of theoretically clean FCNC observables in discriminating new-physics models. 

Generally, inclusive rare decay modes of the B meson are theoretically clean observables. 
For instance the decay width T(B — ► X s ~f) is well approximated by the partonic decay rate 
r(fo — ► 57), which can be analysed in renormalization-group-improved perturbation theory. 
Non-perturbative contributions play only a subdominant role and can be calculated in a model- 
independent way by using the heavy-quark expansion. The inclusive B — > X s £ + £~ transition, 
which is starting to be accessible at B factories |J, represents a new source of theoretically 
clean observables, complementary to the B — > X s 7 rate. In particular, kinematic observables 
such as the invariant dilepton mass spectrum and the forward-backward (FB) asymmetry in 
B — > X s £ + £~ provide additional clean information on short-distance couplings not accessible 
in B - X al . 

Non-perturbative corrections to B — > X s £ + £~ scaling with l/m| and l/m^ can be calcu- 
lated quite analogously to those entering the B — > X s ^y rate [7|, |8], |S], [TTJ. Here the cc resonances 
represent a more serious problem since, for specific values of the dilepton invariant mass, cc 
states can be produced on shell. However, these resonances can be removed by appropriate 
kinematic cuts in the invariant mass spectrum. In the perturbative window, namely when 
me+e- £ m&/2, theoretical predictions for the invariant mass spectrum are dominated by the 
purely perturbative contributions, and a theoretical precision comparable with the one reached 
in the decay B — > X^ is possible. The B factories will soon provide statistics and resolution 
needed for the measurements of B — > X s £ + £~ kinematic distributions. Precise theoretical es- 
timates of the SM expectations are therefore needed in order to perform new significant tests 
of flavour physics. 

In this paper we complete one important step of this program: we present a new calculation 
of QCD-bremsstrahlung and the corresponding virtual corrections in B — > X s £ + £~. This 
effect has already been evaluated in Refs. |Ll], [12| for the dilepton invariant-mass spectrum. 
Here we perform the calculation using a different technique, namely using a full dimensional 
regularization of infrared divergences (both soft and collinear ones). We anticipate that our 
results agree with those of [jy], 12 1 for the dilepton invariant-mass distribution. Moreover, our 
technique allows us obtain the first computation of the soft-gluon (and corresponding virtual) 
corrections in the FB asymmetry: the missing ingredient needed for a systematic evaluation 
of this observable beyond the lowest non-trivial order. Our phenomenological analysis shows 
that these new contributions are rather important, for instance the position of the zero of the 
asymmetry receives corrections of 0(10%). 

The paper is organized as follows: in Section 2 we set up the framework of the calculation; 
we briefly review the status of the various ingredients needed for a systematic analysis of 
the partonic process b — > s(g)£ + £~ in perturbation theory. In Section 3 we present the basic 
expressions needed to describe the spectrum and the FB asymmetry of B — > X s £ + £~ in 
the so-called next-to-next-to-leading logarithmic (NNLL) approximation. The details of the 
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calculation, including a definition of the four-particle phase space and a discussion on the 
regularization scheme, which avoids any ambiguity due to the definition 75, are presented in 
Section 4. Section 5 contains a brief phenomenological analysis of our results, mainly focused 
on the determination of the zero of the FB asymmetry, and Section 6 a short summary. 



2 Theoretical framework 

Within inclusive B decay modes such as B — > X s 7 or B — > X s £ + £~, short-distance QCD 
corrections lead to a sizeable modification of the pure electroweak contribution, generating 
large logarithms of the form a™(m b ) log m (m fe /M heavy ), where M heavy = 0(M W ) and m < n 
(with n = 0, 1,2, ...). A suitable framework to achieve the necessary resummations of these 
large logs is the construction of an effective low-energy theory with five quarks, obtained 
by integrating out the heavy degrees of freedom. With the help of renormalization-group 
(RG) techniques, one can then resum the series of leading logarithms (LL), next-to-leading 
logarithms (NLL), and so on: 

a n s {m b ) \og n (m b /M) [LL], a n s + \m b ) log" (m b /M) [NLL],... (2.1) 

The effective five-quark low-energy Hamiltonian relevant to the partonic process b — > s£ + £~ 
can be written as 

AG 

H ef f = -^fV*V tb £ dQi, M heavy ) O t (fi) , (2.2) 

where 

1 = (s lfl T a P L c) (cYT a P L b) , 2 = (s lfl P L c) {crfP L b) , 

3 = {si»PLb)Y, q {qi*q) , O a = (s / y tl T a P L b) Y, q {(ll^T a q) , 

5 = (s^ lp P L b)J2 q (qrrYq), O e = (s^ lulp T a P L b)E q (qrrYT a q)i 

(2.3) 

Or = -^—m^iscr^P^F^, O s = r|^m 6 (/i) (sa^T a P R b) G% , 

l07T Z \bll Z P 

2 2 

^ 9 = i7r~2 ( s ^&) {i-fi) , d w = {- sl ,p L b) {lYid) 

define the complete set of relevant dimension-6 operators and Cj(/x, M heavy ) the corresponding 
Wilson coefficients. 

Within this framework, QCD corrections are twofold: corrections related to the Wilson 
coefficients, and those related to the matrix elements of the various operators, both evaluated 
at the low-energy scale /1 ~ m b . As the heavy fields are integrated out, the top-quark-, W-, and 
Z-mass dependence is contained in the initial conditions of the Wilson coefficients, determined 
by a matching procedure between full and effective theory at the high scale (Step 1). By means 
of RG equations, the Cj(//, M heavy ) are then evolved at the low scale (Step 2). Finally, the QCD 
corrections to the matrix elements of the operators are evaluated at the low scale (Step 3). 

Compared with the effective Hamiltonian relevant to b —>■ 57, Eq. (|2.3| ) contains the 0(a em ) 
additional operators C 9 and O w . Moreover, it turns out that the first large logarithm of the 
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form \og(mb/Mw) arises already without gluons, because the operator C 2 mixes into Og at 
one loop. This possibility, which has no equivalent in the b — > 37 case, leads to the following 
ordering of contributions to the decay amplitude 

[a em log(m b /M)] <(m 6 ) log n (m fe /M) [LL] , 

[« em log(m 6 /M)] a n s +l {m b ) \og n {m b /M) [NLL] , . . . (2.4) 

Technically, to perform the resummation, it is convenient to transform these series into the 
standard form (|2.1| ). This can be achieved by redefining magnetic, chromomagnetic and lepton- 
pair operators as follows [O, H|: 



16/r 2 ~ q 2 ~ 

O i= _0 4 , <? 4 = -p-Q , (i = 7,. ..,10). (2.5) 

This redefinition enables one to proceed in the standard way, or as in b — > 37, in the three 
calculational steps discussed above |13|, [14]]. At the high scale, the Wilson coefficients can be 
computed at a given order in perturbation theory and expanded in powers of a s : 

Q = CP +* C P + 7 ?Lc? ! > + ... (2.6) 



(4tt) 1 (4 



7T 



Obviously, the Wilson coefficients of the new operators CV-10 at the high scale start at order 
a s only. Then the anomalous- dimension matrix has the canonical expansion in a s and starts 
with a term proportional to a s : 

7 = %(o) + ^L 7 (d + ^L 7 (2) + ... (2 . 7) 

1 4n J (4vr) 2 ^ (4tt) 3/ 1 ' 



In particular, after the reshufflings in ( |2.5| ) the one-loop mixing of the operator O2 with Og 
appears formally at order a s . 

The last of the three steps, however, requires some care: among the operators with a 
non- vanishing tree- level matrix element, only Og has a non- vanishing coefficient at the LL 
level. Therefore, at this level, only the tree-level matrix element of this operator {(Og)) has 
to be included. At NLL accuracy the QCD one-loop contributions to (Og), the tree-level 
contributions to (0 7 ) and (O w ), and the electroweak one-loop matrix elements of the four- 
quark operators have to be calculated. Finally, at NNLL precision, one should in principle 
take into account the QCD two-loop corrections to (Og), the QCD one-loop corrections to 
(0 7 ) and (C10), and the QCD corrections to the electroweak one-loop matrix elements of the 
four-quark operators. 

Let us briefly review the present status of these perturbative contributions to decay rate 
and FB asymmetry of B — > X s £ + £~: the complete NLL contributions to the decay amplitude 



can be found in fl3| , [L4| . Since the LL contribution to the rate turns out to be numerically 
rather small, NLL terms represent an 0(1) correction to this observable. On the other hand, 
since a non- vanishing FB asymmetry is generated by the interference of vector (~ ^7,9) and 
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axial-vector (~ O\o) leptonic currents, the LL amplitude leads to a vanishing result and NLL 
terms represent the lowest non-trivial contribution to this observable. 

For these reasons, a computation of NNLL terms in B — > X s £ + £~ is needed if we aim at the 



same numerical accuracy as achieved by the NLL analysis of B — > X s ^ [15]. Large parts of the 
latter can be taken over and used in the NNLL calculation of B — > X s £ + £~ . However, this is 
not the full story. In particular, the bremsstrahlung corrections presented in this paper for the 
first time are a crucial ingredient, necessary for a complete evaluation of the FB asymmetry 
to NNLL precision. And in the case of the differential rate, the full NNLL enterprise is still 
not complete: 

[Step 1] The full computation of initial conditions to NNLL precision has been presented 
in Ref. [nj. The authors did the two- loop matching for all the operators relevant to 



b — > s£ + £ (including a confirmation of the b — > 37 NLL matching results of |T7 



The inclusion of this NNLL contribution removes the large matching scale uncertainty 
(around 16%) of the NLL calculation of the b — > s£ + £~ decay rate. 

[Step 2] Thanks to the reshufffings of the LL series, most of the NNLL contributions to 
the anomalous-dimension matrix can be derived from the NLL analysis of b — > 57. In 
particular the three-loop mixing of the four-quark operators C?i_6 into O7 and 0% can 
be taken over from Ref. |15|. The only missing piece for a full NNLL analysis of the 
b — > s£ + £~ decay rate is the three- loop mixing of the four-quark operators into Og. In 



Tq] an estimate was made, which suggests that the numerical influence of these missing 
NNLL contributions on the branching ratio of b — > s£ + £~ is small. Interestingly, since 
the FB asymmetry has no contributions proportional to |((9g)| 2 , this missing term is not 
needed for a NNLL analysis of this observable. 

[Step 3] Within the B — > X s 7 calculation at NLL, the two-loop matrix elements of the 
four-quark operator O2 for an on-shell photon were calculated in and quite recently 



confirmed in PUj. An independent numerical check of these results has been performed 



and will be presented in f21| . This calculation was extended in JTT| to the case of an 



off-shell photon (for small squared dilepton mass), which corresponds to a NNLL con- 
tribution relevant to the decay B — > X s £ + £~ . In the dilepton spectrum, this calculation 
reduces the error corresponding to the uncertainty of the low-scale dependence from 
±13% down to ±6.5%. This step also includes the bremsstrahlung and virtual contri- 
butions that are discussed in the present paper. 

In principle, a complete NNLL calculation of the B — > X s £ + £~ rate would require also 
the calculation of the renormalization-group-invariant two-loop matrix element of the op- 
erator Og. Similarly to the missing piece of the anomalous- dimension matrix, also this 
(scale-independent) contribution does not enter the FB asymmetry at NNLL accuracy. 
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3 Basic expressions 

We normalize all the observables by the semileptonic decay rate in order to reduce the uncer- 
tainties due to bottom quark mass and CKM angles: 



r[fe -> X c ev e 



G\m\ 

1927T 3 



V A \ 2 f(z)K{z) . 



Here z = m 2 /ml (m^c denote pole quark masses), f(z) = 
phase-space factor and 

2a s {m b ) h(z] 



1 



k(z) = 1 



37T f(z) 



(3.1) 

\z + 8z 3 — z 4 — 12z 2 In z is the 

(3.2) 



takes into account QCD corrections (the function h(z) has been given analytically in [[22| and 
is quoted in the appendix). The normalized dilepton invariant mass spectrum is then defined 
as 



R(s) 



fT(B^XJ+£-) 



(3.3) 



T(B -> X c eu) ' 

where s = (p£+ + p^-) 2 /m 2 . The other important observable is the forward-backward lepton 
asymmetry: 



A FB (s) 



1 



T(B -> X c ev) J-i 



a cos Op 



ds d cos Of 



-sgn(cos ^) 



(3.4) 



where Qp is the angle between £ + and I? momenta in the dilepton centre-of-mass frame. It was 
shown in || that A FB (s) is identical to the energy asymmetry introduced in [f23fl . 

Following closely — but not exactly — the notation of Refs. [|11], I2| , we present here some 
useful formulae that allow us to systematically take into account corrections beyond the NLL 
level for the partonic contributions to these two observables: 



R(s) 



a. 



4vr 2 



nvts 



V, 



cb 



+(l + 2s' 



4 1 + 



\C, 



u-w^|2 



1 + — T 77 (s 
7T 



\cr(s) 



+ \c 



10 



7T 



-T99U 



+ 12 Re [C 7 new (s)C 9 new ( S )*] 1 + ^r 79 (s) + -5 R (s 



Oi, 



71 



7T 



(3.5) 



A FB (s) 



3a 2 m 


VZVu 


2 


a-*) 2 r 


4tt 2 


v cb 




f(z)n(z) I 



{ S ReK w ( S )*C 9 ncw (^)] (l + ^r 910 (s] 



a 



+2Re[C 1 n ew ( S )*C 7 new (s)] l + ^r 710 (5) +-5 FB (s) . (3.6) 



7T 



7T 



With respect to Refs. (TT|, 0, we have introduced a new set of effective coefficients, defined 
as 



C^ ew (s) 



1 



7T 



<x 7 (s) ) C* 



47T 



s ' s V s ) 



c[ 0) Fi 7 \s) + Ci 0) F?\s) + <5f ^ 
5 



a, , ^ ^ eff , , a. 



C^is) = [l + -a 9 (s))Cf(s)-^ CrFr , (s) + C^Fr(s) + Cr w Ft , (s) 

7T J 47T L 



f ) Fl 9 \s) + Ci 0) Fi 9 \s) + Cf {0) F^i 
CT{s) = (l + %,(*)) 3? . (3.7) 



7T 



The C^ ew have the advantage of encoding all dominant matrix-element corrections, which 
leads to an explicit s dependence in all of them. 

As we shall illustrate in detail in the next section, the terms cxi(s) take into account 
universal 0(a s ) bremsstrahlung and the corresponding infrared (IR) virtual corrections pro- 
portional to the tree-level matrix elements of O7- 10 • The remaining (finite) non-universal 
bremsstrahlung and IR virtual corrections are encoded in rate and FB asymmetry through 
Tj(s) and Sr^b(s). Analytic and numerical results for the Tj(s) and the <7j(s) will be presented 
later on. Here we simply note that the <Jj(s) are defined in order to take into account the truly 
soft component of the radiation, which diverges at the s — > 1 boundary of the phase space. 
With such definition of the (7i(s), the additional finite terms Tj(s) are rather small all over the 
phase space and particularly for large values of s (|rj(s)| < 0.5 for s > 0.3). Substantially 
smaller are the bremsstrahlung corrections not related to Oj_iq <8> O7-10; which are encoded 



in #r,fb("S)- A complete evaluation of 5r(s) can be found in |]12[ , where this effect is shown to 
be at the 0(1%) level. 

The other explicit 0(a s ) terms in ( |3.7|) are due to virtual corrections that are infrared- 
safe. In particular, the two-loop functions F± ^ and the one-loop functions F$'^ correspond 
to virtual corrections to 0\,% and 0$, respectively. These functions have been computed in 
Ref. [11 for small s. The coefficients C^_ lQ) including the one-loop matrix-element contribu- 



tions of 01-6 are defined in close analogy with Ref. |L4| and are reported in the appendix as a 



function of the true Wilson coefficients Cj. Finally, explicit expressions for the latter, evolved 
down at the low-energy scale, can be found in |16| . 



By means of expressions ( |3.5| ) and ( |3.6| ), we can more easily discuss the organization of 
the perturbative expansion for the two observables. According to the arguments in Section 2, 
the formally leading terms are obtained by setting C^ ew = C^ w = 0, neglecting Ti and <5j, 
and identifying Cg cw with the leading term of Cg ff [formally 0(l/a s )}. At this level Ap^(s) is 
clearly vanishing. At the NNL level, when Afb(s) receives its first non- vanishing contribution, 
we should retain the interference of the 0(l/a s ) term in Cg ew with the leading 0(1) terms 
in Cyfo, as well as the sub leading corrections in [Og ew | 2 . Within this approach, the missing 
ingredients for a NNLL analysis of Af-q(s) are only cr 9 and T910. 

As anticipated, the standard LL expansion is numerically not well justified, since the 
formally-leading 0(1/ a s ) term in Og ew is much closer to an 0(1) term. For this reason, it has 
been proposed in Ref. [[□]] to use a different counting rule, where the 0(1/ a s ) term of Og ew 
is treated as 0(1). We also believe that this approach, although it cannot be consistently 
extended at higher orders, is well justified at the present status of the calculation. Within 
this approach, the three 4 new and the two observables [R(s) and v4 FB (s)] are all homogeneous 
quantities, starting with an 0(1) term. Then all <Tj, Tj and 5i functions are required for a 
next-to-leading order analysis of both R(s) and Afb(s). 
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Figure 1: Virtual (up) and real (down) QCD corrections generating the terms and cr, in 
Eqs. (|3.5|) - (|3.7|) . The boxes denote the insertion of £ e fj in (|4.1|) . 

4 Details of the calculation 



In contrast to previous calculations of B — > X s 7 and — ► X s £ + £ matrix elements, we set 
m s = and, following the approach of Ref. [p4| , we employ dimensional regularization to 



regulate all infrared singularities, including collinear divergences arising from the vanishing 
strange-quark mass. 

The diagrams we need to compute are shown in Fig. [l|, where the boxes denote the insertion 
of the following effective non-local Lagrangian 



eff 



Kp 



Cfsi„P L b hn + qfis^PLb £^£ + —Cfsa, u P R bF^ 



wff- 



(4.1) 



where 



Kp 



7T 



V2 



(4.2) 



and, for simplicity, we have omitted to explicitly show the s dependence of Cg S . Using 



instead of the local Hamiltonian in (|2.2|) , we can easily take into account the (two- loop) IR- 
divergent contributions of four-quark operators, whose one-loop b —>■ s£ + £~ matrix element is 
encoded in the non-local coefficients Cf s Q] ( see appendix). 



4.1 Virtual corrections 

The amplitude for the three-body process b(p) — > s(p s )£ + (pi)£~ (p 2 ) can be written as 

2k f r . 



A 



[3] 



m b 



{sP R [m b ^A(s) + (p" + p v s )B(s) + 0(<f )] bI llx Z 
+sP R [m b YA'(s) + (p" + P :)B'(s) + 0(q»)} bl llxl d} 



(4.3) 



where q = pi + P2 and, at the tree level, 



off 



1/^eff 
S 7 



4'(o) _ 



eff 
10 5 



fl'<°> = 



(4.4) 



Employing the same notation, and working in d = 4 — 2e dimensions, the renormalized virtual 
one-loop corrections in Fig. □ are given by 



B 0) 



3tt V a ; Utt, 



eff 



1 5 T- M 



1 + 2s 



2s 



+-Cf 



1 



5 



— - - — - 2Li 2 (s) - 3 ln(l - s) - 10 - 8 In — 



2e 



ln(l - 

+ 0(e) 



3tt v a ; Utt, 



+-Cf 



1 5 

- + - + 2Li 2 (s) 



2e 



Q tt ^ln(l- S ) 

+ ln(l-s) + 10 + 81nf— 



+ 0(e 



a 



^r(i + e)(i- s )- 2 M-^ 



3tt 



a 



4tt 



Ceff 
10 



-3 --Ml 



0(e) 



-^r(i + e)(i- s )- ie 



3tt 



4-7T 



C" 



eff 
10 



1 5 T- ^ 

— o Li 2 is) 

2e 2 4e V ; 



-ln(l-.) 



ln(l 



(4.5) 



The above results include also self-energy corrections of the external legs, not explicitly shown 
in Fig. [|. Effectively, these contributions are included using on-shell renormalization condi- 
tions. For later comparison, we state here the corresponding on-shell wave-function renormal- 
ization constants for the external quark fields, where we explicitly separated ultraviolet (UV) 
and IR poles: 



Z^{m b ) 
Z^(m s = 0) 



1 - 
1 - 



a s 4 / m 2 \ e / 1 



4tt 3 V V 2 
a s 4 / 1 



+ 4 



CUV ClR 



4vr3 Veuv em 



(4.6) 



Because of the conservation of vector and axial currents, UV divergences cancel out completely 
in the terms proportional to C| ff and Cf^. On the other hand, UV divergences proportional to 
Cj S have been eliminated by the on-shell mass counterterm Z mb and by the MS renormalization 
of the corresponding operator using the counterterm Z^: 

a q 16 1 



^77 



1 + 



1 - 



4rr 3 euv 
Og_4 ( m\ 
4tt3 I /i 2 



+ 4 



(4.7) 



S 



After the UV renormalization has been performed, we suppress the /i dependence correspond- 
ing to the IR divergences in order to simplify the notation, as can be seen in Eq. ( |4.5| ). 

Taking into account the full <i-dimensional three-body phase space, virtual IR corrections 



to the rate and forward-backward asymmetry are obtained by means of Eqs. (|4.3| )- (|4.5| ) 
isolating the 0(a s ) terms in 

d d ~ l q 



, r mb d d 1 p s 

[3] 2 (2ix) d - l 2E s (2vr)^ 1 2g 



(2vr) c( - 1 5 d (p -p s -q)x 



x 



d pi d p 2 



1 ^ 1 



spins colours 



(4- 



4.2 Bremsstrahlung corrections to R(s) 

The matrix element of the four-body process b(p) — > s(p s )g(k)£ + (pi)£~ (P2), squared, can in 
general be decomposed as 



spins colors 

where L$ U A denote the leptonic tensors 

1% = tr {fcfkf) and 1% = tr fcfMls) , 
and in addition to s we have introduced the following kinematical variables 



(4.9) 



(4.10) 



J 1, 



k ■ p 



m 



2 ' 



k-p s 



V ■ (pi - Pa) 



2 ' 



k- (pi- p 2 ) 



mi 



mi 



(4.11) 

* 6 m% m% m% 

The asymmetric leptonic tensor vanishes if averaged over the full leptonic phase space 
(symmetric in p\ «-»• p 2 ). On the other hand, in the case of Lg V we can write 

d d ~ 1 pi f d d ~ 1 p 2 



1 s 



(2vr 



27r) d - 1 2£ 1 J (27i) d ~ 1 2E 2 
(l-6)(g 2 )- 



5 d (q- Pl -p 2 )L% 



^23- 4 T(f-e) 



7T 



(4.12) 



Since J5 depends only on q, the variables y and x„ do not appear in the calculation of the 
(symmetric) dilepton spectrum. 

The explicit calculation of the real-emission diagrams in Fig. IT] leads to 



— 2 (?V - qY v ) H 
mi v ' ' 



s 

flU 



IQtt 



K F a s < (1 - s) 



1 l-s 2 M 
^ + + 1-e 



2X r 



P iX.' p iX,' g L * J S 



X a (1 S ) Xr 



1 + 2s - 2es 



°9 



+ 



CefT 
10 



2\ 4 + 2s-4e ~ eS 2 



Cf\ + (6 - 4e)Re (cf Cf * 



+ 



x s (l + 2s) -2x, 



X .11 



°9 



+ 



2^ 4x s (s - 2 + 8x p ) - 8x p (s + 4) 



SXr 



a 



■cff 



4Xo + 8X. 



.1 r 



Re (Cf Cf*) + 0( 



(4.13) 



Infrared singularities arise only by terms proportional to 1/x 2 or l/x s : for this reason only the 
0(e) pieces proportional to these terms (between square brackets) have been explicitly shown. 
Taking into account the full d- dimensional four-body phase space, the radiative differential 
rate can be written as 



dT 



ds 



2 

C s 



(2tt 



,(2-1 



1 



d d ~ l p s 



d d ~ x k 



(2^) d ~ 1 2E s J (27r) d - 1 2E g J (2Tr) d - 1 2q 



d d ~'q 



. 5 i {p _ Ps _ k _ q)I ^ H s 



\3-4e 



where 



16tt 2 

x[l-x^(l-s)]- 2+2£ 

a. 



dx(l - xf-^x 1 - 2 ' / du(l - uu)~ e u 
Jo 



1 — e)m 



5-6e 



7r 2_3e 2 9-10 er (1 _ e)r (|_ e ) r (|_ e ) • 

The integration variables in Eq. (|4.14 ) are defined by 



UJ 



:i - cos^ 



sk i 



X 



2p ■ k 
ml — q 2 



(4.14) 
(4.15) 

(4.16) 



where 8 s k is the angle between gluon and strange-quark momenta in the B rest frame, so that 



x(l — s) 



x s (s,x,uj) 



xuj{1 — x)(l — s) 2 



2 ' ' 2[1 - xu}{\ - ,s)] ' 
Performing explicitly the integrals in x and u we obtain 
dT[4] a.. 



(4.17) 



ds 

x 



x 



4 + - - 4Li 2 (s) - 41n(s) ln(l - s) 
e 2 e 

1 + 2s - 2es 



2 35 1 s(2-3s) , , 



°9 



+ 



U 10 



+ 



4 + 2s - 4e 



{1-8 

+ (6 - 4e)Re ((5f <5, 



eff* 
9 



eff 



+ 



s(4s 2 + 10s - 4) ln(s) - 10s 3 + 7s 2 + 6s - 3 
4(1 -s) 2 

s(6s 2 + 24s - 24) ln(s) - 19s 3 + 18s 2 + 15s - 14 ~ ff 2 

3s(l-s) 2 7 
s(2s + 12) ln(s) - 9s 2 + 4s + 5 



C 



Re (CfC; 



eff* 
9 



(4.18) 



In order to construct an IR-safe observable, dT^/ds has to be combined with the 0(a s ) 
virtual corrections in the non-radiative process. The differential decay rate for the latter can 
be written as 

dT[3] = ^C s s-%1 - S ) 2 - 2 T(1 - e) 



ds 



: {2A 2 [1 + 2s(l - e)] + 4AB(1 - s) + 2B 2 (1 - s) 2 + (A, B A', B')) (4.19) 

10 



in terms of the reduced amplitudes (A, B) and (A', B') defined in ([4.3[). Using Eqs. fl4.4|)-(p~5|) 
and isolating the 0(a s ) terms we then obtain 

J = ^(i-.rrd + «-«) 



X 



1 + 2s - 2es 



a 



cff 



+ 



c 



2 5 , , 2 + 4s, , 
- + - + 4Li 2 (s) + ln(l - s) + 12 



4 + 2s - 4e 



a 



off 



-(6-4e)Re(cfCf* 
ln(l - s 



e s 

4 + - + 4Li 2 (s) + 4 ln(l - s) + 20 + 16 In (JL 

2 5 1 4- 4s 

- + - + 4Li 2 (s) + —— ln(l - s) + 16 + 8 In 



(JL 



(1-8) 



°9 



+ 



°10 



4(4 - s) |cf| 2 - 2(1 + 2s)Re (cfCf*) 



(4.20) 

It is straightforward to check that all divergences cancel in the sum of Eq. ([4.18 ) and Eq. ( |4.20| ) . 
Combining these two equations we can finally obtain explicit expression for the t»(s) and o~j(s) 
functions defined in Eq. ( |3.5| ). 

As mentioned before, we define the universal functions cr 7i9 (s) such that the non-universal 
terms r 77 (s) and t 99 (s) vanish for s — ► 1. This is because in the s — > 1 limit only the soft 
component of the radiation survives and, according to Low's theorem, the latter gives rise to 
a correction proportional to the tree-level matrix. We then define 



09 (s) 



o s 



<t(s) + =■- -In I 



6 



a(s) + - , a 7 (s) 

4 2 2 2 

— Li 2 (s) - - ln(s) ln(l - s) - -tt 2 - ln(l - s) - -(1 - s) ln(l - s) . (4.21) 



3 v ' 3 
With this choice the T;(s) are given by 

2 



9 



9 



?79 (S 



9(2 + s) 
4 

'9(1 + 2s) 



x2, /, n 6s(2-2s-s 2 ), . , 11 -7s -10s 2 
2(1 - s) 2 ln(l - s) + — ^- - — ln(s) + 



(1-s) 2 

2(l-s) 2 ln(l-s) + 3g(1 + s)(1 - 2g) ln(s) 



4(1 -s) 2 4s(3 - 2s) , 
■ ln(l - s) - — r^- ln(s / 



9s 



9(1 



;i-s)2 

^ 2(5 - 3s) 
~ 9(1 -s) 



(1-8) 

3(1 -3s 2 ) 
1 - s 



(4.22) 



and all vanish at s = 1. 

These results are in complete agreement with those recently obtained by Asatrian et al. |TT 
Indeed the u>i(s) of Ref. |Ll] can be written as u> 7 = <j 7 + r 77 /2, u 9 = a 9 + r 99 /2, u 79 - 
(cr 7 + 09 + r 79 )/2. 



4.3 Virtual and bremsstrahlung corrections to Ayb(s) 

A use of naive dimensional regularization is not allowed in the case of the forward-backward 
asymmetry because of the ambiguities arising for d 7^ 4 in the definition of 75 (or in the 



11 



trace of L^)\ in the case of the decay rate, the problematic 75 contribution vanishes because 



of the pi <-> P2 permutation symmetry of the leptonic phase space [see Eq. ( |4.12| )1. To 
circumvent this problem, we employ the following regularization scheme: we treat the Dirac 
algebra of IR-divergent pieces strictly in four dimensions (dimensional reduction). This of 
course simplifies the calculation of the real emission, where only IR divergences appear, but 
it slightly complicates the virtual corrections, where UV divergences also occur. The latter, 
which do not involve any 75 ambiguity, are still computed in naive dimensional regularization 
and renormalized as discussed at the beginning of this section. We stress that, at this level of 
the perturbative expansion, this hybrid regularization scheme is gauge invariant. 

Within the virtual corrections we therefore need to strictly separate IR and UV divergences. 
As a result, the on-shell wave-function-renormalization factors given in Eqs. (|4.6| ) within naive 
dimensional regularization are partially changed. While the IR part of Z^(m&) is not modified, 



0) is different within our hybrid dimensional scheme: 



Z^,(m s = 0) 



1 - 



0^4 / fj^ 

47T 3 \m 
a s 4 / 1 



2c 



47r3 Veuv 



1 2 
+ — 

1 

em 



(4.23) 



Being of pure UV nature, the other Z factors remain unchanged. On the amplitude level, the 
change of Z^(m s = 0) turns out to be the only change of virtual corrections. Taking this effect 
into account, Eqs. ( |4.5|) are modified as follows: 



A 



(i) 

DA 



B 



(i) 

DA 



AW + ^ra + eXi-*)-* ^ 

5 M + |ir(l + e)(l- a )-* fef 



4° 9 + Ys Cr . 



— c c 7 

2s 7 



eff 



A 



'(i) 

DA 



^) + ^r(i + c )(i- a )-* ^ 



4 9 . 



B 



'(1) 
DA 



B >(i) 



(4.24) 



As a consistency check, we have explicitly verified that the finite corrections to R(s), computed 
in the hybrid scheme, are exactly the same as in Eqs. ([4.21 )— (4.22) . To this purpose, the 
bremsstrahlung term with four-dimensional Dirac algebra is obtained by neglecting all the 



0(e) terms in Eq. Q4.13 ). 

The asymmetric part of the squared matrix element of the four-body process 
to Afb(s) — computed with four-dimensional Dirac algebra, is given by 



relevant 



m 



128vr 



+s 



1 

+ - 
s 



-K F CY S < y 



1 1 

— K + 



2x r 



x: 



2yx p - x v (l - s) 
2x p x s (l - s) 

2yx p -x v (l- s) 
2x p x s (l - s) 



XpX s 



(2,(1 



2 

— + — 

ft ft 



2x v ) +y-+^ 2 



x 



Re 



Cw * ( sC f + 2C: 



Re (CfrCf) 



s)) + —{x s + 2s-2xp 



X r 



12 



X y 



2x1 



^1 *^lOC. qXj 'or 



1 — OS )X„ 



Re (c^*Cf 



(4.25) 



As can be noted, H^ U L^ is linear in the leptonic asymmetric variables y and x v . The latter 
can be written as 



— cos 9, 
p ■ k 





2 2 " 


fp-q" 






) 


\ m 2 b , 





1/2 



m 



i , / ,(p-p s )cos# sfc + (p-fc) 
{p-q) q 2 mf l 



(4.26) 



where is the angle between £ + and S momenta in the dilepton centre-of-mass frame and (pe 
the corresponding azimuthal angle. The integrals of y and x v over the leptonic phase space, 
weighted by sgn^os^), can be expressed in terms of s and the two integration variables x 
and u) defined in Eq. ( |4.16|) : 



d d 1 pi f d l p 2 5 (q — pi — P2, 



2E 1 

y(s,x,u) 

x v (s, x, uj) 



2E 2 (2tt, 
(1 - s)P(s,x,u) 
2\l-x(l - s)oj] 



d-2 



sgn(cos#£ 



y 

X i , 



7r l-£ 2 4-4 er (! _ e ) 1 X^(S,X,CJ) 



:i + s - 2x s 



1 - s + 2x s )(l - 2cj) + 4u;x ? 
(l-s) 2 -4(l + s)x s + 4x 2 



(4.27) 



where x SjP are given as in Eq. (|4.17|) and 

P(s, x, uj) = yjl - xuj [4(1 -x) + 2x(l - s)(l - 2a;) + x 2 cj(1 - s)(4 - x(l - Tjj] . (4.28) 

Taking into account also the phase-space integration over the hadronic variables, the regular- 
ized bremsstrahlung contribution to the FB asymmetry can be written as 



dV A 

ds 



where 



d cos 8f 



d 2 T 



[i] 



ds dcos9 t 



-sgn(cos#£ 



C 



16vr 2 



1 S -M _ S )3"4e 



l-2e„l-2e 



dx(l — x) e x 



- u;)~ e ur e [l - xu;(l - s)] 



-2+2e 



m 



ttA t iiv 



y^y(s,x,L)),x v ^x v (s,x,w) 



cr 



m 



5-6e 



7r |-3e 2 io-io e ( 1 _ e ) T fi - e ) r (1 - e) r 



(4.29) 



(4.30) 



Because of the square root in (3(s,x,uj), the integrals on x and uj appearing in (|4.29 ) are 
rather involved and we have not been able to perform them analytically. However, divergent 
contributions arise only in the limit y ^ — (1 — s)/2 from the terms in the first line of ( |4.25| ). 
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Figure 2: The functions fr(s) and fg(s). 



Computing analytically only the latter, we can write 
4a, 



aL [4] 

ds 



3tt 



K 2 F C A S~*(l - Sf-^ 



x 



25 



- + - - 4Li 2 (s) - 4 ln(s) ln(l - s) - tt 2 + — + 



1 3(6 - 7s) 



xRe 



ln(s) 



2 1-s (1-s) 2 
35? (*3f + 2Cf )] +/ 9 (s)Re (sC$*Cf) + f 7 (s)Re (2Cj**Cf) \ . ( 1 .31 ) 



where /7,g(s) denotes regular functions (for s ^ 0), which we evaluate by means of numerical 
methods. The numerical results obtained for /^(s) are shown in Fig. |[ 

The cancellation of IR divergences is obtained by combining ofT^/ds with the C?(a s ) terms 



m 



al [3] 



4CV- e (l - S ) 2 ~ 2 T(1 - e) ^ [16iWj 



ds -, vw ■ (4.32) 

Notice that, as explicitly indicated, in this case one should use the one-loop expressions of A 
and A' computed in the hybrid scheme, given in Eqs. (|4.24j ). Isolating the 0(a s ) terms we get 

^rj^] 1 4a 



is 3/4^(1 - s )-T(i + e) r ( i- e 

2 5 2 -I- 4s 

_ r + ^ + 4Li 2 ( S ) + ^^ln(l- S ) + ll 



Re (sC^Cf) 



14 



+ 



2 5 . . 1 + 5s , . , , / a 

- + - + 4Li 2 s + In 1 - s + 15 + 8 In 

e 2 e s \rrib 



Re (2C^*Cf) ] (4.33) 



and by means of Eqs. ( |4.31| ) and ( [4. 33| ) we finally obtain 



T ™W = 4 + 3(Tb)-||^f ln( S )-i(3-7 S + V)ln(l- s ) + i/ 7 ( s ). (4.34) 

5 Phenomenological analysis 

The results for <Ji(s) and Tj(s) functions, which encode bremsstrahlung and virtual IR cor- 
rections to R(s) and Afb(s), and which we have determined in the previous section, are 
summarized in Fig. |3[ As anticipated, in all cases, except for very small values of s, the 
universal contribution of the <7j(s) is largely dominant with respect to the non-universal con- 
tribution of the Tj terms. f\ The natural scale of these corrections is aia(m b )/n = 0(10%), and 
they thus represent a numerically important effect both in the rate and in the FB asymmetry. 

A detailed phenomenological analysis of R{s) and Ayb{s), including all NNLL effects, will 
be presented elsewhere, together with an independent calculation of the two-loop matrix- 



element corrections [[H]], which extends existing calculations also to the high-dilepton-mass 
region. Here we shall limit ourselves to analysing how the zero of the FB asymmetry (s ) is 
modified by the inclusion of bremsstrahlung and virtual IR corrections. As is well known, this 
quantity, defined by Afb(so) = 0, is particularly interesting to determine relative sign and 
magnitude of the Wilson coefficients C7 and C$. 

Employing the counting rule of Ref. [|ll|], i.e. treating the formally 0(1/ a s ) term of C| ff 
as 0(1) (see discussion at the end of Section 3), the lowest-order value of so (formally derived 
by the NLL expression of A FB ) is determined by the solution of 

s Cf(s ) + 2Cf = . (5.35) 



Using the values of the effective coefficients at fi = 5 GeV reported in Ref. jll| , this leads to 

= 0.14 ±0.02 , (5.36) 

where the error is determined by the scale dependence (2.5 GeV < < 10 GeV). 

Neglecting the small contribution of 5fb [see Eq. ( |3.6|) 1, the next-to- leading order equation 
for s reads 

S0 O 9 ncw ( S0 ) (l + % 910 (s )) + 2O 7 ncw (* ) (l + ^r 710 ( S o)) = , (5.37) 

which leads to 

s nnll = 152 ± 0.OO8 . (5.38) 

3 To simplify the comparison, in Fig. |^ we have plotted the Ti(s) functions and the two combinations 2<7g(s) 
and <7g(s) + 07 (s) which appear in the physical observables. 
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Figure 3: The functions <7j(s) (dashed lines) and Ti(s) (full lines). At s = 0.1 the Tj(s) functions 
are ordered as follows (from top to bottom): r 9 i , r 79 , t 7W , and r 99 (almost identical to t 77 , 
not explicitely shown). 



In this case the variation of the result induced by the scale dependence is accidentally very 
small (about ±1% for 2.5 GeV < \i < 10 GeV) and we believe it cannot be used as a good 
estimate of higher-order corrections: taking into account the separate scale variation of both 
Cg ew and C 7 ew , and the charm- mass dependence, we estimate a conservative overall error on 
s of about 5%. 

We recall that the new effective Wilson coefficients defined in ( |3.7| ) include not only purely 
UV virtual terms, but also universal bremsstrahlung and corresponding IR virtual corrections 
(parametrized in Ui). The global effect of the bremsstrahlung and virtual IR corrections we 
have evaluated is that of reducing by ~ 11% the large (~ 27%) shift between Sq NLL and Sq LL 
induced by the NNLL UV virtual terms:^ 



S s 
so 



Qt 

— M^o) - crg(so) + TYio(so) - t 9 io(s )] = -0-11 • (5.39) 

7T 



On the other hand, it is clear that the reduction of the error (or the reduction of the scale 
dependence) in Eq. ( p.38[ ) is only due to the inclusion of the UV virtual corrections computed 



4 A similar increase of the position of the zero in the FB asymmetry at the ~ 30% level has been found in 
the exclusive B — > V channel |^5| . In the inclusive channel the effect turns out to be somehow reduced by the 
inclusion of the bremsstrahlung corrections, which are of course absent in the exclusive mode. 
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3-1 




Figure 4: Comparison between NNLL and NLL results for Afb(s) in the low s region. The 
three thick (red) lines are the NNLL predictions for \x = 5 GeV (full), and fi = 2.5 and 10 
GeV (dashed); the dotted (blue) curves are the corresponding NLL results. All curves have 
been obtained for m c /m(, = 0.29. 



in Ref . |TTJ . An illustration of the shift of the central value and the reduced scale dependence 
between NNL and NNLL expressions of Afb(s)> in the low s region, is presented in Fig. [|. 
We note that in this region the nonperturbative 1 jm\ and l/m^. corrections to Ap^(s) are 
very small || ^, [T(| and can safely neglected compared to the uncertainties of the partonic 
calculation. 

Interestingly, at this level of precision, the result obtained using the ordinary LL expansion 
is very similar to the one obtained using the phenomenological counting rule of Ref. [ll|. 
Indeed, at a pure NNLL level, i.e. neglecting all 0(a 
a.) terms of Cn ew , we also find as central value Sn NLL = 0.16. 



terms in C° ew and retaining the 0(a s 



x 



6 Summary 

The inclusive B — > X s £ + £~ transition, which is starting to be accessible at B factories ||, 
represents a new source of theoretically clean observables, complementary to the B — > A s 7 
rate. In particular, kinematic observables such as the invariant-dilepton-mass spectrum and 
the lepton forward-backward asymmetry in B — > X s £ + £~ , provide a clean information on 
short- distance couplings not accessible in B — > A s 7. 
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In the present paper we calculated bremsstrahlung and corresponding virtual corrections 
in B — > X s £ + £~. We used a full dimensional regularization scheme of infrared divergences 
(both soft and collinear ones), which also avoids any ambiguity related to the definition of 75. 

For the dilepton-invariant-mass spectrum, these contributions have already been evaluated 
in Refs. fill, |12|1, using a different technique. Our results are in complete agreement with those. 
We also presented the first computation of the soft-gluon (and corresponding virtual) correc- 
tions in the FB asymmetry, with the help of which we evaluated this observable systematically 
to NNLL precision for the first time. 

The new contributions are rather important and significantly improve the sensitivity of 
the inclusive B — > X s l + l~ decay in testing extensions of the Standard Model in the sector of 
flavour dynamics. In particular, the corrections we have computed shift by about 10% the 
position of the zero of the FB asymmetry. The complete effect of NNLL contributions to this 
interesting observable adds up to a 16% shift compared with the NLL result, with a residual 
error due to higher-order terms reduced at the 5% level. 
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Added Note 

After this work was completed, an independent calculation of NNLL corrections to the FB 
asymmetry appeared [p6fl . We find complete agreement with the results of Ref. and 
we thank Christoph Greub for useful correspondence, which allow us to correct a typo in 
Eqs. (|4.26| ) (which affected the s — > behaviour of the plot in Fig. |2|). 



Appendix: Auxiliary definitions 



• The function h(z) describing next-to-leading order QCD corrections to the semileptonic 
decay [see Eq. ( |3.2| )1 is given by p2| : 



h(z) = -(l-z 2 ) (f - f^ + f z 2 ) +zln(z)(20 + 90z-fz 2 + f z 3 
+z 2 \n 2 (z) (36 + z 2 ) + (1 - z 2 ) (f - f z + f z 2 ) ln(l - z) 
-4 (1 + 30 z 2 + z 4 ) \n(z) ln(l - z) - (1 + 16 z 2 + 2 4 ) (6 Li(z) - tt 2 ) 
-32 z 3 / 2 (l + z) [tt 2 - 4 Li(^) + 4 U.{-yfz) - 2 ln(z) In (f^) 

The effective coefficients C^xq appearing in Eq. Q3.7|) are defined as: 



Cf 



a s {jJi) 3 9 3 9 
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Figure Al: Examples of virtual and bremsstrahlung contributions of the four-quark operators 
.6 that are automatically taken into account by the redefinition of the Wilson coefficients 
in Eq. (TQ|) . 



where 



a 



c: 



eff / 



cff 



c 



cff 
10 



+ 
+ 
+ 
+ 



h(z, s) 



An 
An 
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a s {n) 
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4, . . 8 16 z 
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9 V ; 27 9 s 



2 + 



4^ 



42 - s 



X 



2 arctan 



4z-s 



— 2 7T 



for s < 4 2 . 
for s > A 2 . 



(A.l) 



(A.2) 



Note that specific one- and two-loop and matrix-element contributions of the four-quark oper- 
ators 0i_6 (including the corresponding bremsstrahlung contributions) such as the one shown 
in Fig. [AT] are automatically included by the redefinition of the Wilson coefficients C7, Cg 
and C10 given in (|A.1| ). In fact, using this redefinition, the bremsstrahlung and virtual cor- 
rections that were shown in Fig. p] automatically take these effects into account. The Wilson 
coefficients in (|A.1|), which are needed to NNLL precision, are given in p"6|. 
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